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Abstract. We consider cosmological perturbations around homogeneous and isotropic space- 
times minimally coupled to a scalar field and present a formulation which is designed to 
preserve covariance. We truncate the action at quadratic perturbative order and particular¬ 
ize our analysis to flat compact spatial sections and a field potential given by a mass term, 
although the formalism can be extended to other topologies and potentials. The perturbations 
are described in terms of Mukhanov-Sasaki gauge invariants, linear perturbative constraints, 
and variables canonically conjugate to them. This set is completed into a canonical one for 
the entire system, including the homogeneous degrees of freedom. We find the global Hamilto¬ 
nian constraint of the model, in which the contribution of the homogeneous sector is corrected 
with a term quadratic in the perturbations, that can be identified as the Mukhanov-Sasaki 
Hamiltonian in our formulation. We then adopt a hybrid approach to quantize the model, 
combining a quantum representation of the homogeneous sector with a more standard field 
quantization of the perturbations. Covariance is guaranteed in this approach inasmuch as no 
gauge fixing is adopted. Next, we adopt a Born-Oppenheimer ansatz for physical states and 
show how to obtain a Schrodinger-like equation for the quantum evolution of the perturba¬ 
tions. This evolution is governed by the Mukhanov-Sasaki Hamiltonian, with the dependence 
on the homogeneous geometry evaluated at quantum expectation values, and with a time 
parameter defined also in terms of suitable expectation values on that geometry. Finally, we 
derive effective equations for the dynamics of the Mukhanov-Sasaki gauge invariants, that 
include quantum contributions, but have the same ultraviolet limit as the classical equa¬ 
tions. They provide the master equation to extract predictions about the power spectrum of 
primordial scalar perturbations. 
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1 Introduction 

The treatment of linear perturbations has become a prolific field in modern cosmology. The 
advent of inflationary universe models [1] generated great interest in cosmological pertur¬ 
bations and gave rise to a lot of activity in this area. Beyond its early success, the study 
of perturbed systems in general relativity requires great care because one must deal with 
the gauge freedom inherent to the theory, which affects the description of the perturbations. 
Once this problem was realized, most efforts focused on getting a manifestly gauge-invariant 
formalism. 

Nowadays, cosmology is entering a golden age owing to the recent observational progress, 
which has opened new windows to test the predictions of theoretical models [2]. The latest 
observations are providing us with increasingly accurate data of cosmological phenomena 
and, for the first time, it seems possible for astrophysics to think about finding traces of the 
quantum geometric structure of the early history of the universe [3]. For this reason, the 
ultimate hope of the community of physicists working in the quantization of gravity is to 
develop a quantum theory capable of leading to testable predictions. In order to fully capture 
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the quantum nature of spacetime, this theory must involve simultaneously both the geometry 
and the perturbations, with interplay between them. 

Over the last decades, the theory of cosmological perturbations [4], combined with the 
inflationary paradigm [1], has emerged as the framework which conciliates the theoretical 
models of the early universe with observations, since it provides a good approximation to 
the anisotropies of the cosmic microwave background (CMB) and explains quite satisfactorily 
the formation of structures at large scales [5]. The study of the CMB is a powerful tool for 
understanding the universe in its origins. It supports the approximation that the observed 
region is homogeneous and isotropic in a suitable average (demonstrated under certain the¬ 
oretical assumptions [6]). However, this leads to questions about how the anisotropies and 
cosmological structures formed and developed. 

The pioneer work in the analysis of perturbations around classical Friedmann-Robertson- 
Walker (FRW) cosmologies^ is due to Lifshitz [7], as a first modelization of the universe 
considered at a large scale. Nonetheless, it was relatively soon noticed that this analysis had 
been carried out with a specific gauge choice, and hence it did not address the gauge freedom 
satisfactorily. This gauge dependence in the description of the perturbations has caused many 
controversies because keeping track of the gauge modes can get cumbersome, something which 
makes difficult the extraction of the physically meaningful degrees of freedom. An attempt to 
provide a covariant treatment of the perturbations was made by Hawking [8], but this work 
did not resolve totally the gauge ambiguities. It was completed later by Olson [9] for the 
case of an isentropic perfect fluid in a spatially flat spacetime. However, it was Bardeen who 
first constructed a truly gauge-invariant formalism (originally for a perfect fluid), which mixes 
the perturbations of the matter with the perturbations of the four-dimensional metric [10]. 
This work was followed by many other contributions [11, 12]. Likewise, Mukhanov, based on 
Sasaki’s investigations [12], proposed some gauge-invariant held-like variables for the case of 
a scalar held on a spatially hat FRW background, directly related to the co-moving curvature 
perturbations [13]. Mukhanov expanded the action for the gravitational and scalar helds up 
to second order in the perturbations and introduced a gauge-invariant held that completely 
characterized those perturbations and allowed one to rewrite their action exclusively in terms 
of it when the background equations were employed. 

In a standard analysis of primordial huctuations, one studies the perturbations within 
the scheme of Quantum Field Theory (QFT) in a classical and hxed curved spacetime. From 
this viewpoint, the nature of these primordial perturbations is in fact quantum and it is 
rather general to assume that they started out at a very early time as huctuations with a 
small amplitude and gradually grew in time as consequence of gravitational instabilities. For 
doing this, essentially, one represents the perturbations by quantum helds and considers that 
they were initially propagating on a given geometry, such as a de Sitter spacetime, which 
describes rather well the inhationary stage of the universe [14]. Actually, we already know 
from the CMB data that huctuations with amplitudes of 10“^ are sufficient to reproduce 
the cosmic structures observed today. Despite how well this treatment is able to explain 
the present observations, the challenge for quantum cosmology is to build a formalism which 
includes the homogeneous background and the inhomogeneous perturbations and proves to 
be potentially predictive, in order to elucidate whether the relics of the quantum huctuations 
of the early universe may encode information about the quantum character of the spacetime 
geometry itself. 

^These cosmologies are also called Friedmann-Lemaitre-Robertson-Walker cosmologies by many authors. 
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A canonical approach to cosmological perturbations in which the background and the 
inhomogeneities (both in the geometry and in a matter scalar held) were treated together 
quantum mechanically, although hxing the gauge, was developed by Halliwell and Hawking 
[15], in the context of closed FRW universes. They expanded the action up to second order 
in the perturbations and then built the corresponding Hamiltonian. Shortly after, Shirai and 
Wada [16] reformulated this formalism, isolating the quantum dependence on gauge invari¬ 
ants. Actually, they introduced not only a canonical transformation for the perturbations but 
modihed the variables of the background by terms quadratic in those perturbations. Even 
so, physical states still depended on perturbative variables that were not gauge invariants 
(although in a way that was completely determined), and the introduced canonical transfor¬ 
mations were only dehned in a semiclassical limit. Therefore, they depended on the solution 
to the effective equations that was considered in each case. Later, Langlois [17] tried to refine 
and clarify this procedure by using techniques inspired by Hamilton-Jacobi theory in order to 
obtain the gauge-invariant perturbations. Nevertheless, he did not include transformations of 
the background and, besides, obviated the explicit form in which the Hamiltonian constraint 
depends on non gauge-invariant variables, because this should be physically irrelevant. More 
recently, Pinto-Neto and collaborators [18] proposed another approach by means of canoni¬ 
cal transformations which involve the perturbations and the background and, in particular, 
include the Mukhanov-Sasaki (MS) variables. More specifically, they considered the case of 
a perfect fluid and of a scalar field, and reformulated the system so that the global Hamilto¬ 
nian constraint depends only on the gauge invariants and the new background variables. In 
this reformulation, the equations of motion of the background were not used, and the lapse 
function was redefined. Nonetheless, second-class constraints appeared in this reformulation 
procedure which had to be eliminated by reduction and the introduction of Dirac brackets, a 
step that obscures the full gauge invariance of the ultimate system and the role in it of the 
perturbative constraints. 

In the last years, the canonical quantization of cosmological perturbations has received 
as well a lot of attention within the framework of Loop Quantum Cosmology (LQC). LQC [19] 
addresses the quantization of cosmological systems using the ideas and techniques of Loop 
Quantum Gravity (LQG) [20], a non-perturbative and background independent program for 
the quantization of general relativity that provides one of the most promising approaches for 
a quantum theory of the gravitational interaction. Initially, LQC was applied successfully to 
homogeneous models in cosmology, leading to a consistent quantization of the FRW spacetime 
in which the Big Bang singularity is replaced with a quantum bounce, namely the Big Bounce 
[21, 22]. Clearly, the limitation of homogeneity is a restriction that must be surpassed when 
one is interested in studying more realistic scenarios. Therefore, it is natural to try and go 
beyond this simplification and incorporate small fluctuations in the quantum treatment of 
the geometry and the matter content. 

The loop quantization of the FRW background has been combined with a Fock quan¬ 
tization of the perturbations in a so-called hybrid approach in which the whole system that 
describes the perturbed cosmologies is quantized with canonical methods [23~25]. In partic¬ 
ular, in this formalism the perturbed spacetime geometry can be treated as a fully quantum 
entity. This hybrid approach was originally developed in the linearly polarized Gowdy models 
[26]. In this case, the inhomogeneities can be dealt with exactly, rather as pertubatively. An 
alternative proposal, put forward by Ashtekar et al. [27], investigates the propagation of the 
perturbations in the dressed FRW geometry obtained by incorporating the quantum effects of 
LQC on the background. Consequently, in this analysis, known as the dressed metric formal- 
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ism, one loses a full quantum description of the geometry in the perturbed system. Both, in 
the hybrid and the dressed metric approaches in LQC, the perturbations have been expressed 
in terms of gauge invariants, though just after eliminating degrees of freedom in a reduction 
that casts doubt on the covariance of the system. Finally, a third route to analyze the effects 
of LQG in cosmological perturbations consists in assuming corrections in the quantum con¬ 
straints arising from the use of holonomies and the regularization of the inverse of the volume 
operator [28, 29]. Demanding the closure of the modihed algebra of constraints, one can de¬ 
duce the form of those quantum corrections and the corresponding effective field equations for 
the propagation of the perturbations. This strategy is intrinsically covariant, although it still 
rests on a series of hypotheses about, e.g., the form of the possible loop quantum corrections, 
the use of local expansions for non-local quantities, or the independence of the results on the 
existence of superselection sectors. 

In this work, we will provide a canonical formalism for perturbed flat FRW spacetimes 
with a matter scalar held which is specially designed to preserve covariance at the considered 
perturbative level. In particular, the perturbations will be described by the MS gauge invari¬ 
ants, by combinations of the perturbative constraints, and by variables canonically conjugate 
to them. For completeness and self-consistency in the presentation, we will detail the change 
to this description from a classical formulation similar to that of Halliwell and Hawking [15] 
for the case of hat spatial topology. This change will be obtained by means of a canonical 
transformation which will be completed in order to include as well the background variables. 
In addition, we will discuss how to quantize the resulting model using generalized hybrid 
techniques, namely, combining any given quantization of the background with a quantum 
theory for the perturbations (assuming the compatibility of both quantizations in the whole 
system). In this framework, we will show how to extract quantum held equations for the 
MS variables without any gauge hxing. Furthermore, this will be done without appealing 
to semiclassical approximations, like in Refs. [15, 16], nor a Bohm-De Broglie scheme, like 
in Ref. [18]. Instead, we will employ a kind of Born-Oppenheimer ansatz and discuss the 
validity of its application. In this derivation, we will not make use of any background equation 
of motion or constraint, neither at a classical, effective, nor quantum level. Finally, we will 
particularize our discussion to the specihc case of a hybrid quantization in loop cosmology. 

In this sense, our analysis extends other treatments like that of Ref. [25], where the gauge 
freedom associated with the perturbative constraints was hxed, while we do not make any 
such classical gauge hxing here. For that, we will need to abelianize the algebra of constraints 
in our perturbed FRW system. Let us hnally mention that, in the context of LQG as well, 
there have been some attempts to develop a manifestly gauge-invariant perturbation theory 
in the canonical framework by constructing approximate complete observables [30], adopting 
an approach different from ours. 

The plan of the rest of the paper is as follows. Sec. 2 contains the notation and 
the description of the classical perturbed system. In Sec. 3 we will introduce the set of 
transformations for the inhomogeneous modes that leads to gauge invariants, perturbative 
constraints, and their conjugate variables. This will be extended in Sec. 4 to a transformation 
for the whole perturbed system that preserves the canonical structure. In doing this, we will 
have to include quadratic perturbative corrections to the background variables. In addition, 
in this section we will study how these changes modify the zero-mode of the Hamiltonian 
constraint. In Sec. 5 we will explain how to proceed to a (generalized) hybrid quantization 
of the total system formed by the background and the perturbations. Adopting then a 
Born-Oppenheimer ansatz, we will deduce a quantum equation for the propagation of the 


4 


MS variables, discussing its range of validity. Besides, we will derive an effective equation 
governing the evolution of the corresponding classical variables. This formalism will be applied 
to LQC in Sec. 6. We will conclude and summarize in Sec. 7. Finally, we will include two 
appendices with extra details of the calculations. 

2 The perturbed flat FRW model 

In this section we review the classical description of the FRW model with a minimally coupled 
scalar field and with inhomogeneous perturbations (see e.g. Ref. [15]). We consider the case 
of flat spatial topology which (in order to avoid technical problems in the ultraviolet regime, 
as we will comment later on) we assume to be compact, namely that of a three-torus. This 
flat case describes the physically relevant situation in cosmology, since it is in agreement with 
observations, as long as the compactification scales of the three-torus are large compared to 
the radius of our observable universe. The matter content of the model is provided by a 
scalar field $ minimally coupled to the geometry and subject to a potential term. Although 
the discussion can be carried out for generic potentials, we will particularize our analysis 
to the case of a mass contribution for simplicity, explaining just how the main formulas of 
our study generalize for any other potential of the scalar field. On the other hand, as it 
is well known, scalar, vector, and tensor perturbations decouple from each other at leading 
order in the perturbative treatment, and may be considered independently. Actually, for 
scalar fields, vector perturbations are pure gauge and therefore do not contain any physical 
degree of freedom. On the contrary, tensor perturbations are in fact gauge invariant, and 
hence their treatment is relatively simple, at least with respect to the issues of covariance 
of the perturbative formulation that we want to address. Besides, from an observational 
point of view, scalar perturbations are the most interesting ones, since they are ultimately 
responsible for the anisotropies measured in the temperature of the CMB. We will thus focus 
our discussion on scalar perturbations. This model was also considered in Refs. [24, 25] but, 
unlike in this work, the system was reduced by partially fixing the gauge freedom in those 
references. Here, we follow conventions and notation similar to the ones of those works, and 
we refer the reader to them for further details. 

Adopting a canonical 3+1 decomposition, we parameterize the four-metric in terms of 
the three-metric hij induced on the sections of constant time t, a lapse function N, and a shift 
vector (or co-vector Ni). Spatial indices i,j run from 1 to 3. For the unperturbed FRW 
spacetime, these metric functions are completely characterized by a homogeneous lapse A^o(^); 
and by the scale factor of the spatial metric, its square multiplying a static auxiliary three- 
metric ^hij. Instead of using the (positive) scale factor we will consider its (real) logarithm, 
denoted by a{t). On the other hand, we take ^hij as the standard Euclidean metric on the 
three-torus T^, appropriate for the considered case of a compact flat universe. We denote the 
corresponding angular coordinates by 6i, such that 2TT9ijlQ G 5"^, so that the period is Iq in 
each of the orthonormal directions. 

Any function defined on the spatial sections, such as those describing the inhomogeneous 
perturbations, can be expanded in the basis formed by the eigenmodes of the Laplace-Beltrami 
operator compatible with the metric ^hij. In this way, we transform the study of the spatial 
dependence into a spectral analysis in terms of those modes, whose dynamics decouple (in 
fixed backgrounds) at leading perturbative order. Then, as in Ref. [25], we adopt a basis of 
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real Fourier modes, formed by the sine and cosine functions 


Qn,+ {^) 


Qn-{G) 



( 2 . 2 ) 


Here, n - 6 = and n = (u,i,n 2 ,n 3 ) G is any tuple of integers such that its first 

non-vanishing component is strictly positive (this restriction is introduced to avoid repetition 
of modes). These scalar modes have a norm equal to the square root of the auxiliary volume 
Iq of the three-torus, and a Laplace-Beltrami eigenvalue equal to —• u/Iq- 

In the expansion of the inhomogeneities, the vanishing tuple n is not included. The 
zero-modes account for the background homogeneous geometry, parameterized by A^o(i) a-nd 
a{t), and for the homogeneous part of the matter field ‘h, that we denote by (p{t). These 
degrees of freedom are treated exactly at the perturbative order of our truncations in the 
description of the system.^ Using this Fourier expansion, the spacetime metric and the scalar 
field can be expressed as 


hij{t,e) = °hij{9) 


1 + 2 an,e(^)Qn,e(^) 




n,e 

1 lo 


OJi 


r{Qn,e)\ijiG) + - hij{6)QH,e{0) 


N{t, 9) = a 


lVo(i) + e3“W J^5n,.(t)Qn,.(0) 


n,e 


Ni{t,9) = 


= 


a 




+(i) + '^ fn,^{t)Qn,e0) 


(2.3a) 

(2.3b) 

(2.3c) 

(2.3d) 


In these equations, we have defined the constant (T^ = 47rG/(3ZQ), G is the Newton constant, 
the vertical bar stands for the covariant derivative with respect to the auxiliary metric ^hij, 
and e = +,— for cosine and sine modes, respectively. Besides, we have scaled for convenience 
the shift vector and the inhomogeneous part of the lapse function by a power of the scale 
factor^ (and of the mode frequency in the case of the shift). The time-dependence of the 
geometry and matter scalar perturbations is parameterized by the functions 


(2-4) 

In what follows we will omit the time-dependence in these quantities to simplify the notation. 

Following an approach parallel to that in Ref. [15], we can now substitute these expres¬ 
sions in the Hamiltonian form of the gravitational action with a minimally coupled scalar field 

^See e.g. Refs. [15, 16, 18, 31], and the discussion on this issue in Refs. [25, 29[, as well as the full treatment 
beyond perturbation theory adopted in the inhomogeneous Gowdy cosmologies [32[, that confirms that when 
the inhomogeneities are regarded as perturbations, this procedure to deal with the zero-modes is correct. 
®We note that and gn,e(t) are not exactly those of Ref. [25[, owing to the commented scaling. 
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and truncate the result at quadratic order in the inhomogeneous perturbations. In this way, 
one obtains the total Hamiltonian H of the perturbed system at this order of approxima¬ 
tion. As expected, this Hamiltonian is given by a linear combination of constraints, reflecting 
the freedom inherited from general relativity to perform spatial diffeomorphisms and time 
reparameterizations. Specifically, we get [15, 24] 


H = N, 


Hu 


+E«i‘ 


+ E w 


.H' 






1 • 


(2.5) 


Here, denotes the scalar or Hamiltonian constraint of the unperturbed flat FRW model, 
that generates homogeneous time reparameterizations in that system: 

i/|0 = ^(-vr2+vr2+e6“mV). (2-6) 

The constant rh is the mass m of the scalar field conveniently re-expressed as m = ma. 
Besides, we employ the notation TTq to denote the momentum conjugate to any variable q. 
Notice that, in the perturbed case under study, the zero-mode of the Hamiltonian constraint 
(or global Hamiltonian constraint) gets contributions from the inhomogeneities. At our trun¬ 
cation order, these contributions are quadratic in the perturbations. We have included them 
in the sum of terms H^^. For each mode of the perturbations, the contribution is^ 


= - 


-3a 




+ T 


2 / 1 „2 


4,6 + + 4(f 4,6 + 64,6) 

® ^3^n®n,6 “I" 3^nPn,e “I" 3^n®n,6^n,6 “ ^nfn,e^ 

34 (44 “ ^ 4 , 6 ) + ^^an,efn,e + 4,6 }• 


+ e®“m2 


(2.7) 


On the other hand, H^‘^ and are linear in the inhomogeneous perturbations. ^^rises 
from the perturbation of the Hamiltonian constraint around the FRW model in full general 
relativity, while comes from the perturbation of the momentum (or diffeomorphisms) 

constraint. These linear perturbative constraints are given by 


1 r 


|i 21 


H?f = ^ - 27r„7ra._^ 27r,^7r/._^ - (4 + 34)«n,6 - |4e'^“(«n,6 + 

+ + 2/ri,6)], 


1 


H'y = o [ - + ^a(«n,6 + 46fJ,e) + . 


( 2 . 8 ) 

(2.9) 


We note that and do not parameterize physical degrees of freedom, but they are 
instead the Lagrange multipliers associated with these linear perturbative constraints. Finally, 
it is worth remarking that, at the order of truncation adopted in the action, the perturbed 
system is symplectic, with canonical variables given by the zero-modes a and </?, the Fourier 
coefficients = {ofi 6) 6)/n 6} (with I = 1,2,3), and their corresponding momenta 

[15, 24]. 

'^This formula corrects a misprint in Ref. [24], that was not relevant for the results discussed in that work. 
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3 Perturbations in terms of gauge-invariant variables 


As explained in the introduction, our first goal is to describe our perturbative system in 
terms of gauge-invariant variables without fixing any gauge freedom, preserving covariance at 
the level of the perturbative description. With this aim, we first focus our attention on the 
inhomogeneous sector of the phase space which contains the degrees of freedom of the pertur¬ 
bations. In the previous section, these degrees of freedom were parameterized by the variables 
{Xq^’^}, together with their canonically conjugate momenta vr;,^ ttj, ^}, with 

I running from 1 to 3. We will introduce a canonical transformation to describe these per¬ 
turbations in terms of the MS variables, two suitable combinations of the linear perturbative 
constraints, and appropriate conjugate pairs of all of them. Since the intention is to respect 
the canonical structure of the set of elementary variables used for the perturbations, it is clear 
that we need to find, in particular, combinations of the linear perturbative constraints which 
commute, therefore abelianizing the perturbative constraint algebra. Once this abelianiza- 
tion is introduced, the fact that the MS variables are gauge invariant, and hence commute 
with the perturbative constraints, makes straightforward to find a complete set of compatible 
elementary variables for the inhomogeneous sector. The wanted canonical transformation is 
then attained with a convenient choice of conjugate variables. Later on, in Sec. 4, we will 
complete this transformation into a canonical one in our entire system, that is, considering 
not only the inhomogeneities but including also the homogeneous sector of the phase space, 
parameterized by the canonical variables for the zero-modes {a, tTq, (/9, VTip}. In total, the re¬ 
sulting Hamiltonian will be a linear combination of constraints, that are not only first-class 
(as usual), but furthermore that form an abelian algebra. 

So, since in this section we are only interested in the canonical transformations of the 
inhomogeneous sector of our system and in the symplectic structure induced on it, we will 
consider momentarily that the homogeneous sector is a fixed background. With this purpose, 
and only in that sense, we can ignore for the moment the Poisson brackets of the zero-mode 
variables both between themselves and with the perturbations. Hence, in this section we will 
operate with a Poisson bracket in the inhomogeneous sector defined accordingly as 



dF dG 

dx^rdx;t 


dF dG \ 

dx;t dx^r)' 


(3.1) 


where F and G are functions of the perturbations, and might also depend on the background 
variables. 


3.1 Canonical transformation for the perturbations 

Let us start by introducing the MS variables. In terms of the configuration modes the 

modes of the MS field are given by [13, 24, 25] 


Vfia = e 


fn,t H- -{o-n,e + ^n,e) 


(3.2) 


It is straightforward to check that these modes indeed are gauge invariant, as they Poisson 
commute with the linear perturbative constraints: 


}{P) = 0 = }(p). 


(3.3) 






We would like to complete these modes to a set of compatible elementary variables for the 
description of the perturbations. Since the MS variables are gauge invariant, it is natural to try 
this completion by considering combinations of the linear perturbative constraints, with which 
they commute. Besides, in this way, the information about the perturbative constraints will be 
straightforwardly incorporated in our system in terms of elementary variables. In particular, 
imposing them quantum mechanically will be a direct task. The fact that we want to extract 
two compatible variables from the perturbative constraints means, as we have already pointed 
out, that we have to abelianize the algebra of those constraints (and a posteriori, the entire 
constraint algebra of the perturbed FRW system). With this aim in mind, we notice that 
the only non-vanishing Poisson brackets between them is = e^“if|o. It is then 

easy to abelianize their algebra: We only need to replace with the combination 


Hir = m 


= -vr, 


- 3e^“Pf|oa^,, 






<p/n,e- (3.4) 


Actually, this change amounts to a redefinition of the zero-mode of the lapse function. Indeed, 
in the action of the system and up to the quadratic order in perturbations that we are working 
with, we can rewrite the Hamiltonian (2.5) as 


H = N, 


Hu 


+E"i- 


+ ^ 9n,, 


m 


T kfi 




1 ! 


(3.5) 


with the new zero-mode of the lapse function acquiring contributions (of quadratic order) 
from the inhomogeneities: 

No = No + 3e^“ ^ gn,ean,e- (3-6) 

n,e 

Note that the MS variables Vfi^^ remain invariant with respect to the new set of constraints. 
In summary, the set of variables provides a parameterization of the inho¬ 

mogeneous configuration space (inasmuch as that all the variables are compatible) in terms 
of constraints and gauge invariants, as we wanted. 

We now complete the canonical transformation in the inhomogeneous sector by deter¬ 
mining variables canonically conjugate to our new set. It is straightforward to check that 

= 36^,, (3.7) 


is canonically conjugate to n’l, namely = 1, while it Poisson commutes with 

and . Finding the other canonical variables requires a bit more of work. We skip 
the details of the calculation and encourage the reader to check that one can choose 


<^|1 - 


T ^n,e)) 

TTq, 


(3.8) 


as the variable conjugate to the constraint H^‘^, in the sense that 


}(P) 


1, whereas 


+ — (e®“piV/n,e + 37r2 6jj,, 


(3.9) 
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is a momentum conjugate to the MS variable that is ,,, 7r^_ = 1. 


For convenience, we will assign the role of configuration variables to and 

and view the constraints n'l and as momenta (this will simplify the discussion about 
the imposition of the perturbative constraints a la Dirac in the quantization of the system). 
Furthermore, we will use a compact notation for the new set of canonical variables for the 
perturbations, namely we define 



(3.10) 

(3.11) 


In this way, our previous configuration variables and the new ones are related by 
means of a contact transformation. 

3.2 Redefinition of the MS momenta 

In Sec. 5 we will carry out a Fock quantization of the perturbations, and in particular of 
the MS gauge-invariant field. If one reduces the system classically so that it gets described 
by QFT in a curved background, the resulting MS field verifies the Klein-Gordon equation 
of a massive scalar field (with time-dependent mass) propagating in an ultrastatic spacetime 
with compact spatial topology. Remarkably, a series of studies on the uniqueness of the Fock 
representation for Klein-Gordon fields of this type [33-35] proves that the use of the MS field 
to describe the perturbations allows for a unitary implementation of the quantum dynamics 
of the field, and that any other parameterization for the perturbations that includes a time- 
dependent rescaling of the field (as the one employed in Ref. [27]), prevents the dynamics for 
being unitarily implementable in the quantum theory [36]. The unitary implementation of 
the dynamics is possible in a class of (unitarily) equivalent Fock representations with vacua 
that are invariant under the isometries of the three-torus. Moreover, these results require 
a specific choice of momentum for the MS variable among all the canonical pairs, namely 
the evolution can be implemented unitarily as long as the MS modes verify the equation of 
motion 


(3.12) 



at the considered perturbative order. Here the dot denotes, as usual, derivative with respect 
to the time coordinate t. Since the evolution equations are generated by the Hamiltonian 
of the system and the MS variables commute with the perturbative constraints, it is not 
difficult to convince oneself that the above condition on the choice of momenta is achieved 
only if the zero-mode of the Hamiltonian constraint, that contains quadratic contributions 
of the perturbations, does not include any linear term in the momenta conjugate to the MS 
variables.^ 

We want to adhere to the above parameterization of the MS field, and therefore we 
want to choose its conjugate variable so as to eliminate all terms that are linear in the 
MS field momentum from the Hamiltonian constraint. With that aim, we need to modify the 
momentum mode ^ by taking advantage of the freedom in adding to it a linear contribution 


of the MS configuration variable Thus, we introduce a change of the form 



(3.13) 


5 


Actually, the condition happens to be not only necessary, but also sufficient. 
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where F is a function of the homogeneous variables (a, tTq, y?, vr^,). We now have to determine 
this function by appealing to condition (3.12). Notice that, by construction, the new set 

provides a canonical set together with 

In order to find the explicit expression of F we can proceed in a simple way as follows. 
For a moment, we go to the longitudinal gauge, which is picked out by the pair of conditions 
bn,e = 0 and = Tra<in,e + In this gauge, we match the resulting MS momentum 

with the result that was obtained (precisely in that gauge) in Ref. [25]. This procedure 
turns out to provide a unique and well determined answer, showing the consistency of our 
calculations. Moreover, the function F is actually mode independent (therefore our notation 
for it). Its form is 


F = —e + 3— ). 








VTn 


Thus, the modes of the new MS momentum are 


— €■ 

^n,e 


in ^ ' 


— e 


-2a 


+ 7r„ + 3- 


VT,, 




VTn 


(3.14) 


(3.15) 


In order to simplify the notation, in the following we will use again the symbol ^ to denote 
these redefined momentum modes and Ip”’*' for the corresponding set of new momenta. 


3.3 Inversion of the canonical transformation 


For completeness, we finish this section by giving explicitly the expression of the original 
perturbative variables {X”’*"} = {Xq^'', Xp^''} in terms of the new ones { 

The result is 


Ojr) (= 


_ ^ T/n,e 

^a Fg, ^ Vq^ , 


bn.n — 


fn,e = 


Ti'a- = 

^n,e 


+ 


92 

3 ® ’ 

-ayn,e , yn,e 


-—Vpf + ^e^^VpY + — ( + 3 ^'^ 


VTn 


VTn 


VTn 


VTo 


yn,e 

<?1 


I I„4o, ,2 .^2 lT/ri,e 

^^ip 2 ® ) ^2 2 ’ 


■Kh- = 


+ 


/Iq, /Iq, /lo 

I„4a, ,2T/n,e 


- 2-KinT^a + 3— ) K 


vr,. 


VTn 


91 


vr: 




^6q: ~ 2,^T/'fi,e 


(3.16a) 

(3.16b) 

(3.16c) 

(3.16d) 


(3.16e) 

(3.16f) 


11 






4 Full canonical set in terms of gauge-invariant variables 

Recall that in the previous section we regarded the zero-modes, described by the variables of 
the unperturbed FRW model = {a, (p} and {rCp} = {vTc, (a = 1,2), as corresponding 
to a fixed background. Now, we proceed to complete our canonical transformation in the entire 
system, including these zero-mode variables. In this way we will succeed in describing our 
perturbed system in terms of gauge invariants, without any gauge fixing, while keeping its 
full canonical structure. 


4.1 Canonical transformation of the zero-modes 

The action of our system (truncated at quadratic order in perturbations) in terms of the 
original parameterization of the inhomogeneous sector has the form 


S= dt 


E 


W^Wp 


?+ENry 

l,n,e 


n,e 

Pi 




(4.1) 


possibly up to surface terms which are not relevant for our discussion. Here, H{w'^, 
is the total Hamiltonian expressed in terms of the original variables employed in Sec. 2, 
{u;“} = {Wq,Wp} and = {X^’^,Xp/*^}, as in Eq. (3.5). Let us focus on the Legendre 

term that contains the information about the symplectic structure of the full system, that we 
call W: 


W = 


o. Lft.e 


n,€ vn,e 
Ql Pi 


(4.2) 


Our goal is to find a canonical transformation relating the previous variables with new zero- 
modes w°‘ = {Wq,Wp) such that W retains its canonical form when expressed in terms of the 
gauge invariants for the perturbations, namely 


W = 


^ a l^n,E 


n,eT^n, 
Ql ^Pl 


(4.3) 


Instead of reproducing all the details of the lengthy calculation that allows one to deduce 
the form of the desired transformation (and which is essentially based on the consideration 
of our change of perturbative variables as a canonical transformation in the inhomogeneous 
sector that depends on a series of time-dependent external variables, describing the homoge¬ 
neous degrees of freedom), we will only sketch the main steps in the derivation. First, we 
notice that the relations X”’*^ = given in Eqs. (3.16) do not mix different modes 

and are linear. Therefore, the original variables X”’*^ can be expressed in the following way 


xr = 


E 

m 


n.e 


dVn 


+ 




Qm 


dVp^ 


K 

n,€ Pm 


(4.4) 


where the partial derivatives are functions of the zero-modes w°‘ and the frequency a;„ only 
(and m = 1,2,3). Taking this relation into account in the derivatives with respect to time 
that appear in W, performing several time integrations by parts, disregarding total time 
derivatives which contribute at most by surface terms at initial and final times (assuming 
a well posed variational principle for the system), and truncating up to quadratic order 
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in perturbations, as well as recalling that is a canonical set as long as one ignores 

the homogeneous sector, one can obtain the expression of the new canonical homogeneous 
variables as functions of the old ones and of the new perturbative variables (or 
of if preferred). In this way, one gets new zero-mode variables which are given by 
the original ones plus some additional contributions that are quadratic in the perturbations. 
This result can be inverted at the considered perturbative order, allowing one to express 
the original zero-mode variables as functions of the new phase-space variables for the entire 
system (zero-modes plus inhomogeneities). With this procedure, one arrives at 


l,n,e 

l,n,e 

In these expressions we have to understand the original variables for the perturbations as 
functions of the new ones as given by Eqs. (3.16), with {w°‘} = {a, tTq,, y?, tt,^} replaced 
with {«)“} = {d, tTq,, (^, in those equations. Note that this replacement is consistent 
within our truncations, since old and new zero-modes differ in terms that are quadratic 
in perturbations. In this sense, let us also clarify that the partial derivatives in the above 
identities must be taken keeping constant. In Appendix A we give the explicit expressions 
of the zero-modes w°‘ in terms of the new phase-space variables w°‘ and 


ft y-’’ 

vn,e 


5 , 7 ,. 


dX,r 

5,7,5 


y-n,e 


5,7,5 


(4.5a) 

(4.5b) 


4.2 Hamiltonian in terms of gauge-invariant variables 

Once we know the canonical transformation relating the original phase-space parameteriza¬ 
tion, with the new one that uses gauge invariants for the perturbations, {u)“, 

we can obtain the form of the Hamiltonian in the new formulation. In particular, let us con¬ 
sider the zero-mode of the Hamiltonian constraint. In the original variables, this zero-mode 
is given in Eq. (3.5) by the sum 

+ (4.6) 

n,e 


which contains the contribution of the homogeneous FRW model and the correction quadratic 
in perturbations. 

Since the difference between the original and new variables for the homogeneous sector is 
precisely of quadratic order in the perturbations, the substitution of the relations between the 
old and new variables for our system leads to the following new expression for the zero-mode 
of the Hamiltonian constraint: 


^ dw^ 


b 


(4.7) 


where the difference w°’ — w°‘ must be regarded as a function of the new phase-space variables, 
given by the last terms in the two expressions (4.5a) and (4.5b). Note that this difference is a 
sum of independent contributions of each of the inhomogeneous modes, so that we can write 


n,e 


(4.8) 
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Besides, evaluation of and at the new zero-mode variables w°‘ is done in Eq. (4.7) 
by simply replacing the original phase-space coordinates directly with these new ones in 
the functional dependence of the considered contributions to the constraint. To derive the 
above expression, it suffices to expand H\q in series around the new zero-mode variables and 
recall that we are truncating the total Hamiltonian at quadratic order in the perturbations. 
From these considerations, we see that the quadratic contribution of the perturbations to the 
zero-mode of the Hamiltonian constraint in our new variables takes the form 'Yhfie -^ 12 '*^ > with 


= 


Hf'- + Y: Au-S, 


3g|i 

: Q^a 


(4.9) 


Indeed, this is the form that one would expect if one regarded the change of variables for 
the perturbations as a canonical transformation that depends on a series of time-dependent 
external variables (namely, the zero-modes) with dynamics governed by the Hamiltonian H^q 
[17]. 

If one carries out the calculation explicitly, one gets the following quadratic contribution 
to the constraint of our system: 


77“’" = 77“’" + 7^2 "7^10 + + F'y - 3 


3q: 


9 


TAn,e I \ TAn,e 

* no ' ^ no I 


TTn 


P2 


P3 I '^pa ' 


(4.10) 


In this expression we have defined 


77”’" = 


iOn + e -I- 


1 -b 15(^2 - I2(p'^ - 

VTa vri 






1 - 93 + ) (H”’")^ + - 97r|) {yf) 


TT 


n,e\2 


,-2q 


- 3- 


vr 


(vr? -b 37r|) Ticp + e^°'vr?(p'Ka 


p ^ _ 


1 e 


a L 
—4a 


.VT, 


Yn,epn,e _ o " V’ 


<?l <12 


Tin 


Pi 92 ’ 


2 TTn 


67r^H”’"+(6e3“77|o + 57ri)H^ 


n,€ 

92 


— 4q 


TT? 


/^6q:~2~^ , 0^3 0.^2^ \ , ,,2T/fi,e ^—6a^ 

ye m ipTla “b 37r^ 1tT^TT(pJ Vq^ F OJ^Vq^ ^e 


- 3e 


—2a '^‘P Tm,e 


VT, 


V . 

pi 


(4.11a) 

(4.11b) 

(4.11c) 

(4.11d) 


In the total Hamiltonian, we can account for the term ^ F^Hiq at the considered pertur¬ 
bative order by means of a new redefinition of the zero-mode of the lapse function, with the 
new one given by Nq = iVo(l -b e -^ 12 ’*^) 1®®® ^9- other terms in Eq. (4.10) dif¬ 

ferent from 77^’" contribute to the perturbative linear constraints and their presence amounts 
to a redehnition of the corresponding Lagrange multipliers, that we now denote as and 
77^ g. In summary, we obtain at our quadratic order a total Hamiltonian constraint of the 
form 


77 = IV, 


77i, 


+EF 


+ 


yn,e 
E '^P2 




n,e 
P3 ‘ 


(4.12) 
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In Appendix A we give the explicit expressions of the original lapse function N and shift 
CO- vector Aj in terms of the new phase-space variables and Lagrange multipliers 

As we see, the terms provide finally the contributions quadratic in perturbations to 
the zero-mode of the Hamiltonian constraint in the constructed gauge-invariant formulation. 
As expected, these terms are precisely those that depend exclusively on the MS variables 
and their momenta, and hence they are gauge-invariant quantities. For obvious reasons, we 
call their sum the MS Hamiltonian. Moreover, as we had anticipated, these terms contain 
no linear contribution of the momenta of the MS variables. In fact, the only contribution 
from these momenta is quadratic, and its mode-independent coefficient is constant up to a 
power of the scale factor. We notice as well that the expression given for is linear in the 
momentum of the zero-mode of the scalar field, This linear expression has been obtained 
by taking advantage of the identity vr^ = 2A|oe^“ + ~ that can be used in Eq. 

(4.10) at the coast of redefining the zero-mode of the lapse function at the considered order 
in perturbation theory. Thanks to this, in Sec. 5 we will be able to interpret the zero-mode 
of the Hamiltonian constraint in a certain approximation as a Schrodinger-like equation that 
dictates the quantum evolution of the inhomogeneities in some family of states. 

5 Hybrid quantization and Born-Oppenheimer ansatz 

We can now proceed to quantize the symplectic manifold which describes our cosmological 
system, and to impose the classical constraints a la Dirac, i.e., as operators that annihilate 
physical states in the quantum theory. We recall that, in our classical formulation, we have 
split the phase space in homogeneous and inhomogeneous sectors, starting with the expan¬ 
sion in modes given by the eigenfunctions of the Laplace-Beltrami operator on the spatial 
sections. The homogeneous sector describes the zero-modes, and its degrees of freedom can 
be parameterized by the canonical variables {rc“} = {ci, tTq,, (^, On the other hand, the 
inhomogeneous sector describes all the non-zero modes of the perturbations, and can be pa¬ 
rameterized by that is a set specially selected to make manifest the 

covariance of the formulation at the considered perturbative level. In this way we have pre¬ 
pared the classical system to now carry out a (generalized) hybrid quantization, in which 
one adopts specific quantizations for both the homogeneous and inhomogeneous parts. On 
the tensor product of the corresponding representation spaces, one imposes the constraints 
quantum mechanically, so that the total system is not at all trivial. 

In fact, it seems natural to assume that there exists some regime of the quantum dy¬ 
namics, in between a fully quantum gravity regime and the regime where QFT in hxed curved 
backgrounds makes sense, in which the most relevant quantum effects of the geometry are 
those affecting the zero-modes, and that perturbations, while still having a quantum descrip¬ 
tion, can be represented in more standard ways, like e.g. as in QFT in curved spacetimes. 
Hence, we can adhere to a hybrid quantization based on this hypothesis, adopting a Fock 
quantization for the MS gauge invariants (and possibly for the rest of variables for the per¬ 
turbations) and, in principle, a different type of quantization, of quantum gravity nature, for 
the homogeneous sector. Note that this hybrid approach is different from treating the pertur¬ 
bations as a test held propagating in a quantum background, which is the viewpoint adopted 
in the dressed metric proposal in Ref. [27]. Rather, the inhomogeneities and the homoge¬ 
neous sector, even being represented with different quantization methods, form a symplectic 
manifold that is quantized as a whole. The Hamiltonian constraint of the system affects both 


15 


sectors, encoding the backreaction of the perturbations on the homogeneous background inas¬ 
much as this is kept in the quadratic truncation of the action. 

In this section we will maintain a general discussion without specifying the concrete 
quantization of the homogeneous sector. Moreover, most of the details are easily extensible 
to other quantizations of our MS variables different from the Fock (or QFT-like) quantization 
that we adopt here. We will simply assume that the quantization of the zero-mode variables 
provides a representation for the canonical commutation relations such that the operators for 
the homogeneous FRW geometry commute with those representing the homogeneous sector of 
the matter held, and in turn all of them commute with the elementary operators representing 
the variables for the inhomogeneities. We will denote (8> the kinematical Hilbert 

space for the homogeneous sector, such that the operators representing the homogeneous 
FRW geometry are dehned on while the operators for the homogeneous matter sector 

are dehned on . Finally, concerning the quantization of the homogeneous matter sector, 

we will also assume for simplicity that the operators that represent functions of (p act by mere 
multiplication, while the operator representing Ti(p acts as a generalized derivative. 

In previous applications of this strategy of a hybrid quantization [23-26], the polymeric 
representation of LQC was adopted for the background geometry, together with a standard 
Schrodinger representation for the zero-mode (p of the matter held. Although we will briehy 
explain the particularization of our study to such a representation in the next section, we 
intend to make here a much more general discussion, by no means restricted to the framework 
of LQC. Therefore, our analysis generalizes the quantum treatment presented in Ref. [25] 
in two important ways: a) instead of reducing the system classically, we carry the linear 
perturbative constraints to the quantum theory, with no gauge hxing; and b) we do not 
adopt here any specihc representation for the homogeneous sector, and in particular for the 
FRW geometry. In this sense, the formulation can be adapted to any quantization proposal 
in FRW cosmology. The physics resulting from each specihc proposal might be discussed and 
compared afterwards. 

5.1 Fock representation for the perturbations 

For the inhomogeneous modes describing the MS gauge invariants, we adopt here a Fock 
representation similar to that of Refs. [24, 25]. Remarkably, this quantization is selected 
by the criteria of: a) invariance of the vacuum under the spatial isometries, and b) unitary 
implement ability of the dynamics in the regime in which one recovers a QFT in a curved 
background for the MS held (in any hnite time interval) [33, 34]. As mentioned in Subsec. 
3.2, these criteria remove the ambiguity in a possible rescaling of the MS modes by functions 
of the homogeneous variables, and they assure a unitary implementability of the evolution 
and a standard quantum mechanical interpretation in the regions where a QFT in a (generally 
effective) background is recovered. Actually, these results require the spatial sections to be 
compact; that is why we choose our hat model to have three-torus topology. More concretely, 
the above criteria pick out not just a Fock representation, but a family of them such that 
they are all unitarily equivalent [34]. In particular, this family contains the representation 
in which the annihilation-like variables, Ufj., and the creation-like variables, at , are those 
naturally associated with harmonic oscillators of frequency ojn- A more specihc choice of Fock 
representation among the privileged family selected by those criteria can be done on the basis 
of further requirements not just on operators that are constructed from exponentials of the MS 
variables and their momenta (namely, the Weyl algebra), including linear operators, but on 
other physically relevant operators like e.g. the held Hamiltonian. It seems natural to demand 
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that the MS Hamiltonian be well defined quantum mechanically (and be essentially self- 
adjoint). Additional conditions related with regularization may be important, although one 
would expect that regularization schemes should be a direct consequence of the quantization 
of the system, since this includes now the background, rather than having to import them 
from QFT, where one usually appeals to them claiming the absence of a quantization of the 
curved spacetime. 

Assume then that we take a particular Fock representation for the MS modes in the 
family picked out by the commented criteria of vacuum invariance and unitary dynamics 
in the standard QFT regime. Let us call and dt^ the corresponding annihilation and 

creation operators, acting on the Fock space F, such that [d^^£,dt, ^,] = . A basis of 

states is provided by the occupancy-number states |AA), where AA denotes an array of (positive 
integer) occupancy-numbers, one for each mode. The creation operator acting on these states 
excites the mode labeled by ft and e, increasing the corresponding occupancy-number in one 
unit, as usual. 

5.2 Quantum representation of the constraints 

In our gauge-invariant formulation the classical constraints Poisson commute, and therefore 
they can be imposed quantum mechanically without troubles as far as their quantum counter¬ 
parts commute as well. We assume that our quantization satishes this property and impose 
them independently. 

Let us start with the linear perturbative constraints, which classically are Vp-’^ = 0 for 
I = 2,3. These constraints are straightforward to impose adopting a quantization, for the 
part of the inhomogeneous sector parameterized by (again with I = 2,3), such that 

the operators representing the momenta act as derivatives (or as translations, in integrated 
or discrete versions) with respect to the configuration variables V^’^. Then, the constraints 
amount to require that physical states are independent of those variables. Hence, after im¬ 
posing these constraints, we can just restrict the discussion to a representation space of the 
form Ti = 0 A" to study all physically relevant quantum states. Let us empha¬ 

size that the restriction to these states is quantum mechanical, and not a classical reduction. 
Besides, let us note that these states are not yet fully physical, since we must still impose the 
zero-mode of the Hamiltonian constraint. 

More challenging is the imposition of this global Hamiltonian constraint, for which we 
will only be able to provide approximate solutions. Let us first focus on its quantum repre¬ 
sentation. Classically this constraint is given by 

H^o + Hi2^e-^^H = 0, (5.1) 

with Lf|o defined as in Eq. (2.6) but in terms of the background variables (d, tt^, (^, vr,^), and 
iL |2 = defined in Eq. (4.11a). In what follows we will focus on the densitized 

Hamiltonian constraint® H = 0. 

Before continuing, let us introduce some convenient notation. We first define 

= tt? - (5.2) 

®We could also choose to work with the non-densitized constraint and carry out the densitization at the 
quantum level (see e.g. Ref. [25]). 
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so that the contribution of the homogeneous sector to H is = (vr^ — 'H^^^/2. In 

addition, we introduce the following functions of the zero-modes: 


and we call 


■do = — 

'^OL 

= e2“, 

■dl = e“2“7r| -h 

= ^19 - 18^^ + (l - 2<f^) , 


0n,e ^ + ^9)(yn,.)2 ^ 

n,e n,e 


(5.3a) 

(5.3b) 


(5.3c) 


(5.4a) 

(5.4b) 

(5.4c) 


Then, 2e^".ff|2 = — (©e+Qo^Ti^)- Let us assume that we can represent the quantities introduced 

in Eqs. (5.2) and (5.3) as densely defined operators, '&o, '^e, and -df, on (8) 

(while acting as the identity on J^), and the objects given in Eqs. (5.4) as densely defined op¬ 
erators, ©o and ©e, on Likin'"Then, we obtain the following operator representing 
H: 



- K 


( 2 ) 


©. 



where we have adopted the symmetrization 


(5.5) 


(^0o7r<^^^ = ^ ^©oTT,^ + vf^0o) = ^[fr<^,0o] + (5.6) 

5.3 Born-Oppenheimer ansatz 

We will now analyze solutions to the (zero-mode of the) Hamiltonian constraint of the system 
by adopting the ansatz 

^ = T{a,(:p)'tp{J\f,(f), (5.7) 

where the dependence on the MS variables is denoted by the label Af of the occupancy-number 
states for the inhomogeneous subsector and a denotes the dependence on the geometry 
of the homogeneous FRW sector.^ Thus, the ansatz for the quantum states is that their 
wave functions T can be separated into two factors, one which depends on the homogeneous 
degrees of freedom of the FRW geometry, and the other on the MS gauge-invariant modes. 
The ansatz allows for states that present different rates of variation in these two factors with 
respect to the zero-mode (p of the matter scalar field, and when this happens the Hamiltonian 
constraint can be approximated and simplified. Hence the name of Born-Oppenheimer for 

^This notation does not imply that we adopt a representation for the homogeneous geometry in which 
we have an operator representing & that acts by multiplication, but it is rather a symbolic way to indicate 
functional dependence on the homogeneous geometry sector. 
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the introduced ansatz. From this perspective, it is convenient to interpret as an internal 
time for the system (at least in some intervals of the evolution). We recall that we are 
assuming a representation of the homogeneous matter sector such that the functions of (p act 
as multiplicative operators on 

Let us for a moment consider the unperturbed FRW model. Its Hamiltonian constraint 
operator is proportional to vr? — . In the evolution picture that we are employing, where 

(p plays the role of the time, positive-frequency solutions of this constraint are given by 

r(d, p) = U (d; (p)x{a), (5.8) 


'' ( 2 ) 

U being the corresponding unitary evolution operator. If "Hg is self-adjoint, then one can 

project on its positive part (P.P.) and take y^P.P.(?ig^^) as the (^-dependent self-adjoint oper¬ 
ator that generates the time evolution in Eq. (5.8). To maintain the discussion more general, 
and since the self-adjointness of T-Cq might not be guaranteed, or a straightforward definition 
of a square root for it might not be available, we simply assume that the FRW part of the 
quantum states has the form (5.8), where the family of evolution operators tj is such that 

there exists a self-adjoint operator Hq satisfying [vr,^, U]U~^ = Hq. Note that the concepts of 
self-adjointness and unitarity that we use here are those corresponding to the FRW-geometry 
part of our Hilbert space, * Actually, within our perturbative framework, it is natural 

to consider F as an approximate solution to the Hamiltonian constraint of the homogeneous 
sector, and not necessarily an exact one. In fact, as we will see, F will supply an approximate 
solution of the FRW model inasmuch as the square of T-Lq is close to "Hg in the sense that 

{HoY — ?ig^^ + [vr^,T^o] is negligibly small on F. 

On the other hand, in Eq. (5.8) we normalize the state x(d) to the unit in . This 
part of the quantum state can be understood as the state of the FRW geometry at a fixed 
initial time (pQ (and hence it is independent oi (p). A natural election would be to choose 
a very semiclassical state x, so that it is highly peaked on a certain homogeneous geometry 
(and, if this is possible, that it remains peaked under the evolution dictated by U). 

Let us plug the ansatz (5.7) in the constraint equation® .H'L = 0. Taking into account 

that 


= r(7r,^V’) + (^or)V’5 (5.9a) 

= r(7r|V') + 2(4r)(^^V’) + + {(4)'r}V’, (5.9b) 


the constraint can be rewritten as 


{ {{no? - 77® + [7r^,77o])r}v- + 2(77or)(7r^V’) + - ^o, 0o](rV’) 

- {0e + (0o7^o)s}(r^) - 0o{r(fr^V’)} = 0. (5.10) 


®In a conventional Schrodinger representation in which = —idcp we have the standard unitary evolution 
operator 

U{a-,p) = V 


exp (i f dipilo{a,ip) 


The symbol V denotes time ordering with respect to ip. For simplicity, we have fixed the reduced Planck 
constant h equal to the unit. 

®Since one expects that solutions do not belong to the kinematical Hilbert space, one should rather impose 


the constraint on some kind of generalized states (ipl with the adjoint action, in the form = 0. We do 

not do so because we do not want to complicate the notation even more. 
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We note that,with our assumptions on the representation of the zero-mode of the scalar field, 
the operators [■7r^,^o] a^nd [tTi^,©,,] depend on if, that acts by multiplication as an operator, 
but are independent of -kf. Besides, notice that the first term of the above equation is the 
correction coming from the consideration of wave functions F that are not exact solutions to 
the constraint of the homogeneous sector (to see this, it suffices to make -0 constant and the 
operators ©o and ©e of the MS modes equal to zero). 

Now, let us consider the approximation that, regarding the homogeneous FRW geometry, 
only the terms corresponding to the expectation values on F are relevant. In other words, 
when taking the inner product of the left hand side of Eq. (5.10) with F with respect to the 
FRW geometry (namely in we disregard possible quantum transitions from F to other 

states mediated by the action of the constraint. From Eq. (5.10) it is not difficult to see that, 
for this approximation to hold, we only need the following operators to have small relative 
dispersions on the state F for all values of (p (or, in other words, the following operators must 

be highly peaked on F along the evolution in ip)\ i) T-Lq, ii) 'de, hi) —^dp'&o + {'&oT~io)s + 
and iv) —idpHQ + where we have defined^*^ 


- idpd = -ilo,6], (5-11) 

for any operator O. In principle, the condition on the operator iv) may be satisfied with a 

suitable choice of T-Lq, in agreement with our comments above. 

Assuming that the considered approximation is valid, we then get 


- (0o)r) 


(©e + {QoHo) g)r + i{dfHo — -di^0o)r + — {'Hof) 


V’- 


(5.12) 


Here we have introduced the symbol (0)r to denote the expectation value of the operator O 
on F in the Hilbert space , namely, only with respect to the inner product of the FRW 
geometry. Note that, in general, the result is an operator acting on (8) J~. 

The above equation leads to a generalized Schrodinger equation for the evolution (in p) 
of the inhomogeneities provided that the following conditions are satisfied: 

3-) (0o)r has to be negligible as compared to ('Ho)r- This is valid as long as the quadratic 
contribution of the MS modes given by ©o remains small when compared to the intro¬ 
duced generator of the (^-evolution in the FRW case, something that certainly fits in the 
perturbative scheme that we have adopted for the treatment of the inhomogeneities. 

b) k'^'if; has to be negligible as well. The self-consistency of this assumption is checked in 
Appendix B. 

c) i{dfi-Lo — ^dipQo)r has to be negligible in comparison with (©e+ (0o’Ho)^)r- Otherwise, 
the Schrodinger equation would contain a term destroying the unitary evolution of 
the MS modes and which is not present in the classical field equations of these gauge 
invariants. 

the conventional case = —id(p, d,^ is the total derivative in the Heisenberg picture corresponding to 
a time evolution in ip generated by T-Lq. 
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Indeed, if these three conditions are satisfied, we get the generalized Schrodinger equation 

^ (0e + (©o^io)^)? + - (?io)^)r , 

TT^'ip = - 2 —X- Ip. (5.13) 

2(?io)r 

On top of the previous conditions, let us consider as well the requirement 

d) — ('Ho)^)r has to be negligible in comparison with (©e + (0o’Ho)^)r- This is 

natural to assume (once (d<^'Ho)r has been ignored) if T is an approximate solution of 
the unperturbed FRW model, as argued before. 

If the four conditions a)-d) hold, then the Schrodinger equation simplifies to 


e + (©o^o) 5 

2(?io)r 

Correspondingly, the Hamiltonian generating the dynamics of the perturbations in the internal 
time ip is given by (©e + {&o'Ho) g)r/‘^{'Ho)r- This Hamiltonian is just the MS Hamiltonian, 
with its dependence on the variables of the homogeneous geometry evaluated at the expec¬ 
tation values corresponding to the quantum state F, and divided by the expectation value of 

Ho on it. 

We see that assuming the Born-Oppenheimer ansatz and introducing a controlled num¬ 
ber of approximations allows us to recover standard QFT for the gauge-invariant pertur¬ 
bations, which propagate on an effective homogeneous geometry than can be regarded as 
dressed. The term dressed refers here to the fact that this geometry is not a classical one, but 
it retains the main quantum corrections to the geometry that the MS variables feel [27]. 


-V’- 


(5.14) 



5.4 Effective equations for the Mukhanov-Sasaki variables 

Employing the Born-©ppenheimer ansatz and the approximation that the state F remains 
highly peaked on the operators that encode the effect of the FRW geometry in the zero-mode 
of the Hamiltonian constraint (and hence, strictly speaking, without the need to introduce 
the rest of approximations that lead to a Schrodinger equation for the wave function of the 
MS modes), one can further derive effective classical equations for the MS variables. These 
effective dynamics treat the perturbations as classical, namely they replace the annihilation 
and creation operators of the gauge-invariant inhomogeneities by their classical counterpart. 
The assumption that this replacement is acceptable in order to get effective equations does 
not seem too stringent, because the zero-mode of the Hamiltonian constraint is quadratic 
in the MS gauge invariants. Since the functions of p act by multiplication, we can see the 
quantum dynamics of the MS modes ruled by that Hamiltonian constraint as corresponding 
to time-dependent harmonic oscillators, with p playing the role of internal time. 

In the previous subsection we saw that the evolution of the MS variables is governed by 
Eq. (5.12), which can be interpreted as the result of imposing a constraint of the form: 

Cper = + Drip)TTp + Er{p) - (Oe + (©O^^o)^ “ (5.15) 

Here, Hr and Ey are two functions of p which depend on the state F of the homogeneous 
geometry, and which we do not specify because they are irrelevant for our calculations. As 
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for the momentum tT;^, it is supposed to act as a generalized derivative with respect to ip, 
and we neither need to specify it at this level of discussion. The constraint operator Cper is 
imposed on the sector of the model composed by <8) -T. The corresponding classical 

evolution generator, denoted by Cper, is then obtained by replacing by its explicit form as 
a generalized derivative, and the MS operators and Vp-,’^ by their classical counterparts, 
according to our comments above. 

Taking into account the densitization of the constraint [set in Eq. (5.1)] and the def¬ 
inition of the homogeneous part of the lapse function, one can check that Cper/2 generates 
reparameterizations in a time T that, at leading perturbative order, is related with the proper 
time by dt = ae^°‘dT. Moreover, we can change to a conformal time r/r, adapted to the dressed 
FRW geometry associated with the state T. Observing expressions (5.4), we see that all the 
dependence of Cper/2 on the MS momenta 1^’^ is given by the term (^9e)r(4p”’'^)^/2 coming 
from (0e)r- Then, it is natural to define^^ lodrjr = {'d(,)rdT. This change of time is well de¬ 
fined because (i?e)r is a c-number (depending on (p through T), and it is monotonous inasmuch 
as the operator i/e representing should be positive. It is worth emphasizing that 

we would have failed to dehne a change of time parameter had the time derivative dijr/dT 
been an operator. Hence, the expectation value on T is essential in order to introduce the 
above change of time. We also point out that the change is state dependent, and hence the 
properties of the evolution in the times T and rjr can be quite different when considered in 
the physical Hilbert space of the system. 

In summary, the evolution of in the time ryr is given (under Poisson brackets) by 




2(R 


e/r 


{^r>Cper}, 


(5.16) 


where dp^ denotes the derivative with respect to r/r. The evolution of is simply given by 
dp^Vq^''' = l^Vp^''. Using this result and taking again the time derivative we get the following 
effective MS equations: 


(f 
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+ (^e + ^o)r 


where we have defined = IqUj'^ and 


° (^e)r ■ 


/fl \ _ ,2 ((^o'Ho)s - 

(^o)r = Iq --y--. 

(^e)r 


(5.17) 


(5.18a) 


We note that the last term in the square brackets of this MS equation is a function of only (p, 
and hence of time, when the scalar held is evaluated on the solutions to the effective equations. 
This factor contains quantum modihcations with respect to the standard MS equation, which 
encode the most relevant quantum geometry effects of the homogeneous background. The de¬ 
rived equations are still of harmonic oscillator type with time-dependent frequencies. Besides, 
no dissipation term appears and the equations are hyperbolic in the ultraviolet regime, where 
dominates in the square brackets. We have included the contribution of (^o)p) although, 
in view of our discussion in Appendix B, we expect it to be negligible in practice as it is 
proportional to fh?, in the studied case of a mass term as the potential. 

is possible to see that, with our choice of numeric factors, this definition of conformal time is not 
sensitive to the choice of Zq. 
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Let us say that, although in our analysis we have assumed for simplicity that the potential 
of the scalar field is given by a mass contribution, actually we can easily extend the discussion 
to a general potential W(cp). For that, we simply need to replace in Eq. (5.3c) the two powers 
of with 2W{(p), and the term fh? with W"{(p), whereas the factor m^ip in Eq. (5.3a) has 

to be replaced with W'{(p) (here the prime denotes derivative with respect to 0). Then, for 
the approximations related with the Born-Oppenheimer ansatz to hold, instead of requiring 
a small mass, we need to require small variations of the potential. 


6 LQC representation for the homogeneous sector 


So far in our discussion we have left unspecified the concrete quantization adopted in the 
FRW-geometry part of the homogeneous sector. For the sake of illustrating our analysis with 
an example of particular interest, in this section we will adopt the polymeric quantization 
employed in LQC. There are plenty of references where the details of the polymeric quantiza¬ 
tion of the unperturbed FRW model can be found (see e.g. [19, 21, 22, 37]). More specifically, 
we will adhere to the so-called improved dynamics prescription of LQC [ 22 ] and, among the 
possible symmetric orderings for the Hamiltonian constraint operator, we will adopt the pre¬ 
scription put forward in Ref. [37], as it proves to be most convenient both for theoretical and 
practical purposes [38]. Moreover, in addition, for the homogeneous scalar field, we will use 
the standard Schrodinger quantization which is usually taken in LQC. Besides, in the per¬ 
turbed model, when representing the operators of the homogeneous sector that are coupled 
to perturbations, we will follow Ref. [25], that adopts the factor ordering of Ref. [37] as well. 

In the loop formalism, the geometry is described by the Ashtekar-Barbero su{2) connec¬ 
tion and by the densitized triad [20], that form a canonical pair. In FRW cosmologies, owing 
to homogeneity and isotropy, they are respectively determined simply by two dynamical vari¬ 
ables, c and p, with Poisson bracket equal to SnCj/S, where 7 is the Immirzi parameter [39]. 
The mentioned improved dynamics scheme accounts for the existence of a minimum non¬ 
vanishing eigenvalue A of the area operator in LQG. This scheme involves a transformation 
to the new variables 

|p|3/2 



The sign of p determines the orientation of the triad. The new variables verify { 6 , v} = 2. 
The variable v is proportional to the volume of the homogeneous model, which is finite for 
the three-torus spatial topology under study, and given hy V = 27rG7\/A|u|. These variables 
for the homogeneous geometry are related with those of previous sections, (d, tTq), via the 
canonical transformation 


e = 


/ 37VA 
I 2c7 ' 


1/3 


tTq, = 


( 6 . 2 ) 


On the other hand, for the zero-mode of the matter scalar field, it is convenient to 
introduce the following scaling by a constant, (f) = ^0 = ^ 0 ' since this is 

the usual parameterization employed in the LQC literature. 

In LQC one polymerizes the connection, something which means that the connection 
coefficient c has no well-defined operator in the quantum theory but one represents instead 
its holonomy elements, given by the exponentials of b. As a consequence, the FRW-geometry 
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sector of the kinematical Hilbert space, , is the span of basis states |r;), with r; G R, 

normalized with respect to the discrete inner product {v'\v) = 5^. The operator v , which acts 
by multiplication, v\v) = v\v), has a discrete spectrum. We are taking the reduced Planck 
constant h equal to one (as in Sec. 5.3) to simplify the notation. The holonomy operators 
N±p_ = produce a constant shift in the label of these states, N±fi\v) = |?;± 1 ), 

as one can deduce from the commutation relations [A^^, u] = i{Nfi,v}. As a result of this loop 
representation, the classical expression cp = 2'KGjvb gets promoted in the quantum theory 
to a symmetric version of vsm(b) multiplied by 2ttG^, where sin( 6 ) = i(iV_ 2 g — N2p)/2. We 
choose the symmetric ordering proposed in Ref. [37], after which the operator representing 
cp becomes 



2VA 


sgn(u)sin( 6 ) + sin( 6 )sgn(u) 


yi/2_ 


(6.3) 


As we have said, for the zero-mode of the scalar field, we adopt a standard Schrodinger 
representation, with = id^j, and (j) acting by multiplication, so that = T^(R, dc/)). In 


A /2') 

total, this LQC representation yields the following expression for the operator "Hq ' [see Eqs. 
(5.2) and (5.5)]: 


rt-o — 




dvrG ydvrGq- 


(6.4) 


The operator Oq annihilates the zero-volume state \v = 0) and leaves invariant its orthogonal 
complement. Moreover, it leaves invariant the subspaces Tif formed by states with support on 
the semilattices Cf = {±(e-|-4n)|n G N}, where e G (0,4]. are separable, in contrast with 
the original Notice also that, in each of these sectors, the homogeneous volume v has a 

strictly positive minimum e (or negative maximum —e). Below we will represent the quadratic 
contributions of the inhomogeneities to the (zero-mode of) the Hamiltonian constraint by 
operators that also leave these semilattices invariant, that therefore get superselected. In the 
following, we will restrict the discussion, e.g., to spanned by states with positive v G 

Let us consider now the representation of the homogeneous terms entering the quadratic 
contribution of the inhomogeneities to the zero-mode of the Hamiltonian constraint, namely 
the operators 'do, '&e, ^■iid idg representing the ??-terms in Eqs. (5.3). These terms are affected 
in principle by some quantization ambiguities. As it was done in Refs. [24, 25], we will 
introduce a symmetric factor ordering that tries to respect, as far as possible, the assignations 
of representation made in the homogeneous sector of the system as follows, i) We represent 
the inverse of the volume with the standard regularization of LQC, namely [1/H] is the cube 
of the regularized operator 



A'irG'yy/A 




(6.5) 


This operator annihilates the state \v = 0) and commutes with V, and hence it is well defined 
on the subspaces Hf. ii) The products of the form where / is an arbitrary function, 

are represented with the symmetric factor ordering (/((/>)vr 0 -|- 7 r 0 /((^))/ 2 . iii) We adopt an 
algebraic symmetrization in factors of the form V''g{cp), that are promoted to the operators 
y'’/ 2 ^yf/ 2 , -vvrhere g is any function and r a real number. This algebraic symmetric factor 
ordering is also adopted for powers of the inverse volume, iv) Even powers of Qq = 
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are represented by the same powers of the operator Oq, as in FRW. And v) in the case of 
odd powers of Oq; l^t us say with k equal to an integer, we adopt the representation 

l^ol^^ol^ol^) where |fio| is the square root of the positive operator fig and Aq is defined 
exactly as rio in Eq. (6.3), but with holonomies of double length. In this way, its action 
only shifts v in multiples of four units, and hence preserves the superselection sectors of the 
homogeneous geometry. Following these prescriptions, we obtain the operators 
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(6.6a) 

(6.6b) 


(6.6c) 


which are densely-defined on Tif (8) L^(M, (i(^). These results coincide with those of Ref. 
[25] except for some factors of the inverse of the volume operator, which appear now as the 
inverse powers of the volume operator itself. This difference arises from our different choice 
of densitization for the zero-mode of the Hamiltonian constraint at the classical level (the 
densitization was done at the quantum level in that mentioned, previous work). 

As we have already said, the validity of the Born-Oppenheimer approximation depends 
both on the particular representation chosen for the above homogeneous operators, and on the 
properties of the homogeneous states F. We refer the reader to Ref. [25] for more comments 
at this respect regarding the loop quantization. 


7 Conclusions 

In this work we have developed a covariant formulation of the gravitational system that 
describes a perturbed flat FRW spacetime with a compact three-torus topology and minimally 
coupled to a scalar field. In our perturbative scheme, we truncate the action of the model 
at quadratic order in the perturbations, and focus our attention on scalar perturbations. 
We expand the spatial dependence in Fourier modes, using a basis of eigenfunctions of the 
Laplace-Beltrami operator of the spatial sections. Besides, we treat the degrees of freedom 
of the zero-modes exactly, up to the order of our truncation. The total Hamiltonian of the 
resulting model is a sum of constraints, formed by a linear combination of the zero-mode 
of the scalar (or Hamiltonian) constraint, and of all the inhomogeneous modes of two local 
constraints which are linear in the perturbations. One of them comes from the perturbations 
of the Hamiltonian constraint and the other corresponds to the perturbation of the momentum 
(or diffeomorphisms) constraint. The zero-mode of the Hamiltonian constraint is, in turn, 
formed by the constraint of the unperturbed model plus contributions that are quadratic in 
the perturbations. 

The first important aspect of our treatment is that we do not fix the gauge freedom asso¬ 
ciated with the linear perturbative constraints. In order to deal with them, we abelianize their 
algebra. We achieve this by replacing the linear perturbative Hamiltonian constraint with a 
suitable linear combination of it and of the constraint of the unperturbed model. At the per¬ 
turbative order of our truncation of the action, this replacement amounts to a redefinition of 


25 











the zero-mode of the lapse function. As a result of this process, we are able to parameterize the 
perturbations of the model with these Abelian linear perturbative constraints and the modes 
of the MS held, together with their corresponding canonically conjugate variables. Since the 
MS modes Poisson commute with the constraints, we indeed attain a parameterization for 
the perturbations fully adapted to gauge invariance. 

Another goal of our work consists in extending the canonical transformation from the 
perturbations to the entire system, namely considering not only the inhomogeneities, but 
also the homogeneous sector of the model (the zero-modes describing the degrees of freedom 
present in the unperturbed FRW model). In this way we get at our level of truncation in the 
perturbations a fully covariant description of the entire symplectic manifold, and not just of 
its inhomogeneous sector formed by the perturbations, something that would have required 
treating the zero-modes as fixed, rather than as genuine degrees of freedom. This is an 
important step towards the quantization of the system, inasmuch as without this completion 
of the canonical transformation, the system would have lost its symplectic structure. As 
a result of this canonical transformation, the zero-mode of the lapse gets redehned again, 
absorbing a quadratic contribution of the perturbations, but keeping the original freedom as 
an unspecihed Lagrange multiplier. In addition, also the Lagrange multipliers accompanying 
the perturbative linear constraints get redehned similarly (though in this case with linear 
contributions of the perturbations), so that the hnal Hamiltonian is a linear combination of 
(hrst-class) commuting constraints. We have derived the explicit expression of the spacetime 
metric in terms of our gauge-invariant phase-space variables and new Lagrange multipliers. 

In our construction, we have started with a form of the linear perturbative constraints 
that has been suitably scaled by powers of the scale factor. It is worth explaining that we 
might have started as well without this scaling of the perturbative constraints: The difference 
would have simply led to another dehnition of the momentum of the zero-mode of the FRW 
geometry in the final canonical set of variables for the entire system. 

Once the classical description has been completed, we have proceeded to quantize the 
system. For this, we have proposed a hybrid approach that combines any quantization of 
the homogeneous FRW sector of the model with a more standard field description of the 
inhomogeneous degrees of freedom. The utility of this approach rests on the assumption that 
the most relevant quantum geometry effects are those affecting the zero-modes of the system, 
and therefore may require a more careful analysis, while perturbations, even being quantum, 
can be described in a more conventional way. In other words, the philosophy behind the 
hybrid quantization consists in adopting a representation for the homogeneous sector of the 
model capable to capture its quantum gravity nature, while the perturbations are treated 
essentially along the lines of QFT in curved spacetimes. 

We have followed the Dirac approach of imposing the constraints as quantum operators 
in order to find physical states of the system. As the classical constraints Poisson commute, 
it is natural to pass to quantum counterparts that commute as well, and hence we can 
impose them independently. The linear perturbative constraints generate translations in 
their canonically conjugate variables. Their quantum imposition then means that physical 
states do not depend at all on those gauge degrees of freedom. Therefore, thanks to our 
abelianization of the constraint algebra, we see that the process of imposing the perturbative 
linear constraints in the quantum theory leads to quantum states that depend only on the MS 
modes and the homogeneous sector. Remarkably, this is precisely the kind of states obtained 
with the approach adopted e.g. in Ref. [25], even if the avenue followed in that case included 
a gauge hxing of the perturbations, although the description of the remaining physical degrees 
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of freedom was done only in terms of gauge invariants. In this sense our results, in the present 
work, support those obtained with gauge hxing and show the covariance of the conclusions 
attained previously, at the considered perturbative level. After the imposition of the linear 
perturbative constraints, the only constraint on the model is the zero-mode of the Hamiltonian 
constraint. 

To solve this Hamiltonian constraint, we have adopted an ansatz for physical quantum 
states of Born-Oppenheimer type, and discussed the validity of such approximation. This 
Born-Oppenheimer ansatz has a similar motivation as the one proposed in Ref. [40]. It 
regards the zero-mode of the scalar held as an internal clock and assumes that the MS modes 
do not affect much the motion of the zero-modes of the geometry. The ansatz assumes a 
separation of the quantum dependence on the FRW geometry and on the MS gauge invariants. 
If the state of the FRW geometry remains highly peaked with respect to three operators [the 
operators listed as i)-iii) in the paragraph above Eq. (5.12)] when it changes in the internal 
time, and this evolution hts suitably that of the unperturbed system, the zero-mode of the 
Hamiltonian constraint provides a quantum dynamical equation for the MS variables in terms 
of this time. Furthermore, under certain conditions, the ansatz allows one to introduce a series 
of approximations that lead to a Schrodinger equation for the MS modes, characterized by a 
Hamiltonian that is given by the MS Hamiltonian corrected by the fact that its dependence 
on the homogeneous geometry gets dressed with quantum corrections. In this manner, one 
recovers a QFT for the gauge-invariant perturbations, that propagate in a fixed (and in general 
non-classical) spacetime. 

Generically, when one is dealing with quantum fields in curved spacetimes, there is am¬ 
biguity both in the way that one parameterizes the held (since one can introduce scalings by 
background functions) and in the Fock representation chosen for that held. Remarkably, for 
the system under study, there exist uniqueness theorems that pick out a preferred param¬ 
eterization (in fact, a canonical pair for the held) and a class of unitarily equivalent Fock 
representations for it [33-36]. These choices are selected by the natural criteria of having a 
vacuum invariant under the symmetries of the background and a unitary implementation of 
the dynamics. Our parameterization for the MS modes matches precisely this description.^^ 
No regime with a unitary QFT would have been accessible if we had chosen a different pa¬ 
rameterization for the MS held (by means of a time-depending scaling). 

Another important conclusion of our analysis is the robustness of the class of dynamical 
equations that govern the evolution of the MS gauge invariants. We have shown that the 
dynamics of these invariants are always given by a harmonic oscillator equation, with a 
frequency that depends on the internal time, as long as the Born-Oppenheimer states are 
highly peaked with respect to the operators that feel the FRW geometry in the zero-mode 
of the Hamiltonian constraint, and irrespectively of whether this constraint equation can be 
approximated by one of Schrodinger type. To arrive at this conclusion, essentially our only 
additional assumption has been that the effective dynamics of the MS modes can be derived 
with the substitution of annihilation and creation operators by their counterparts as classical 
variables. The MS equation obtained in this manner contains quantum modifications to the 
time-dependent part of the frequency, but the modification is the same for all modes, and 
in particular it does not affect the ultraviolet behavior of the classical equations for these 
perturbations in general relativity. 

^^Nonetheless, we emphasize that our strategy of a hybrid approach to the quantization can be easily 
adapted to other parameterizations. 
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Let us also note that, in order to obtain the quantum evolution for the perturbations, 
we have not employed a semiclassical approximation for the homogeneous geometry. We do 
not even need to exactly solve the unperturbed model to determine the quantum state of the 
background geometry. Indeed, within our perturbative scheme, it is enough that we consider 
approximate solutions, whose difference with respect to the exact ones can be neglected 
perturbatively. 

We finish by emphasizing that, if we believe in QFT in curved spacetimes, it is natural 
to think that there is a deeper quantum regime, before reaching a complete quantum gravity 
description, in which our hybrid quantization makes sense. This hybrid approach encodes the 
main quantum gravity effects and it is suitable to potentially predict whether there exist mod¬ 
ifications of quantum gravity nature regarding cosmological observables, for instance in the 
power spectrum of the CMB. It therefore offers a framework to extract physical consequences 
of quantum gravity in cosmology. 
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A Spacetime metric and scalar field in the gauge-invariant formulation 
A.l Zero-modes 

In this appendix, we give the explicit expressions of the original canonical variables for the 
homogeneous sector of the phase-space of our system, {re®} = {a, tTq,, <y9, in terms of the 
new canonical set introduced in the text, formed by the variables w°‘ and Substituting 

Eqs. (3.16) and their corresponding derivatives into Eqs. (4.5) yields the desired expressions: 
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A.2 Lapse function and shift vector 

The reformulation of the model in terms of gauge invariants for the perturbations leads to the 
Hamiltonian given in Eq. (4.12), with Lagrange multipliers {A^o, and constraints 

expressed in terms of the new zero-mode variables {d, tTq, vr,^} and new variables for the 
inhomogeneities Comparing Eq. (4.12) with Eq. (3.5), and taking Eq. (4.10) 

into account, we get (within our perturbative truncation scheme) 
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From these expressions we deduce the inverse relations (again within the order of our pertur¬ 
bative truncation): 
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where ^ is given in Eq. (3.16a). These formulas dehne the modes of the lapse function 
and the shift vector, given in Eqs. (2.3b) and (2.3c), in terms of the variables and Lagrange 
multipliers of our gauge-invariant formulation, up to terms that are neglected at our truncation 
order in the perturbations. 

To complete the expressions of the spacetime metric components in our formulation, 
dehned by Eqs. (2.3), as well as of the scalar held (2.3d), we just need to substitute 
and g by their values given in Eq. (3.16), and to replace a and in terms of a and (p [see 
Eqs. (A.la) and (A. 2a)]. Recall that the difference between old and new variables for the ho¬ 
mogeneous sector is quadratic in the perturbations, and therefore some contributions arising 
from the above replacement may be disregarded as negligible within our perturbative trun¬ 
cations (namely, if they provide quadratic or higher order corrections to the inhomogeneous 
modes, or cubic or higher corrections to the zero-modes). 
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B Negligible second derivatives: self-consistency of the approximation 


In this appendix we investigate the self-consistency of neglecting the second derivative of the 
wave function of the MS modes with respect to the homogeneous part of the scalar field in 
the constraint equation of the Born-Oppenheimer states. 

Suppose that, in Eq. (5.12), we can ignore (0o)r compared to ('Ho)r- This assumption 
is justihed on the basis of our perturbative scheme. Then, from that equation, we get 
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Now we act on this equation with and compute the value of the second derivative of ?/>. 
Taking into account that [vr^, (0)r] = —i{^(pO)r and eliminating terms which are negligible 
perturbatively, we arrive at 
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Starting with this equation, it is possible to see by iteration (operating repeatedly with 
Tt(p) whether it is legitimate to assume that the action of on ip is negligible compared to 
the action of vr^, and in particular if this happens for n = 1. For this, one can see that it is 

sufhcient that — ('Ho)^)r be negligible compared with linear terms in the perturbations, 

and that (d^^T^Q^^ — d^(?io)^)r be negligible as well compared to terms of quadratic order, all 
this assuming that 

i) (d^0e + d^(0o?io)5)r is negligible with respect to terms of quadratic order in the 
perturbations. 


ii) {d^T-Lo)r and (d<^0o)r are at most of the order of quadratic terms. 

Let us recall that —idipO = [vr^, O] — [T-Lq, O]. The hrst term is a generalized derivative 
with respect to the explicit dependence on (p of the considered operator O. For the relevant 
operators in our discussion, this dependence comes exclusively from the potential of the matter 
scalar held, that we have considered to be given by a mass term. Then, if the possible values 
of the mass are considerably small, the term ([vr^, 0])r may be treated perturbatively, e.g. by 
expressing the mass value as a certain power of the amplitude parameter of the inhomogeneous 
perturbations. 0n the other hand, when O is any of the theta-operators, the second term 
(['Ho)0])r gives a non-vanishing contribution, whose details will depend on the particular 
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representation chosen for the homogeneous geometry, and on the properties of the state F. 
In consequence, whether the above conditions are true or not will depend on these particular 
details and should be carefully checked. 
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